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Abstract 



We present the results of QCD analysis of out-of-event-plane momentum distribution 
in 3-jet e + e~ annihilation events. We consider the all-order resummed perturbative 
prediction and the leading power suppressed non-perturbative corrections to the mean 
value (K out ) and the distribution and explain their non-trivial colour structure. The 
technique we develop aims at improving the accuracy of the theoretical description 
of multi-jet ensembles, in particular in hadron-hadron collisions. 
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1 Introduction 



Physics of e + e~ annihilation into hadrons is being used as a testing ground for developing QCD 
analyses of multihadron production in hard processes. The studies of the structure of hadron jets 
produced in e + e~ annihilation have led to establishing the standards for the accuracy of QCD 
calculations. The state-of-the-art includes sophisticated perturbative (PT) analysis and taking 
into account the leading non-perturbative (NP) effects. The latter show up as contributions 
suppressed as a power of the hardness scale Q. In particular, in the distributions and means of 
various jet shapes the NP effects are typically O (1/Q). 

At the level of the PT description of e + e~ jet-shape observables, these standards include |jj 

• all-order resummation of double- (DL) and single-logarithmic (SL) contributions due to 
soft and collinear gluon radiation effects, 

• two-loop analysis of the basic gluon radiation probability and 

• matching the resummed logarithmic expressions with the exact O (a^) results. 

The NP technology developed in recent years |2| unambiguously provides the exponent of the 
leading power correction to a given observable. To quantify, in a universal way, the magnitude 
of the genuine confinement contribution one has to 

• carefully address the problem of merging PT and NP contributions, and 

• take into consideration two-loop effects in the magnitude of the NP correction, in particular 
those due to non-inclusiveness of jet observables (Milan factor). 

These are the today standards for the QCD predictions concerning typical e + e~ hadron systems, 
that is two-jet events. 

QCD treatment of hard processes involving hadrons in the initial state such as DIS and, 
especially, hadron-hadron collisions, seldom (if at all) reaches such an accuracy. At the same 
time, hadron physics generated by ensembles of more than two hard partons is theoretically quite 
rich. Because of QCD string/drag effects, the structure of accompanying flows of relatively soft 
hadrons depends on the geometry and colour topology of the multi-prong hard "parton antenna" . 
This makes it also practically important: measuring spatial distribution of hadron flows provides 
a tool for identifying the nature of the underlying hard collision on event-by-event basis ||, 
markedly at LHC. 

NP effects should also be sensitive to the geometry and colour structure of the multi-prong 
hard-parton system. This issue, as far as we can tell, has not yet been addressed in the literature. 
As a first step in this direction, we report in this letter the results of the QCD analysis of the 
distribution of three-jet e + e~ annihilation events in out-of-event-plane transverse momentum 
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We consider the near-to-planar "3-jet region" 

T~T, / »T m = ^, (1.1) 

where T,Tm, and T m are the thrust, thrust major and thrust minor, respectively. The thrust 
axis and the thrust-major axis are set equal to the z- and y-axis respectively. The {yzj-plane is 
defined as "the event plane" . 

At parton level the events in the region ( |1 . 1| ) can be treated as being generated by a system of 
energetic quark, antiquark and a gluon accompanied by an ensemble of secondary (soft) partons. 

We study the distribution of events in the cumulative out-of-plane transverse momentum 

3 n 

Rom = y^ bosi + \kjx\, (1.2) 

a=l i=l 

with pi,p2, P3 the momenta of the hard partons generating three jets, and fc» the secondary parton 
momenta. 

Similar to 2-jet configurations in the region of two narrow jets (1 — T, Mjj/Q 2 , C, B>t,w <S 1), 
physics of near-to-planar 3-jet events, in the region T rn <C 1, involves DL (Sudakov) multiple 
radiation effects. 

Soft gluon radiation dominance makes it possible to resum DL terms to all orders. The PT 
answer for the -fT ut-distribution is expressed in terms of (the exponent of) the "radiators" — the 
basic probabilities of a soft gluon emission off each of the three event-generating hard partons. 
Subleading SL corrections due to hard collinear radiation in one of the jets can be embodied into 
the hard scale of the corresponding radiator. SL effects due to soft inter-jet particle production 
induce the dependence on the event geometry (angles between jets), both in the PT and NP 
contributions. 

At Born level there are only three massless partons with momenta P a (Pi > P2 > P3), so 
that K OVLt = 0. The most energetic Pi lies along the thrust axis, and we define the second 
most energetic momentum P 2 to have a positive y-component. Kinematics (energies and relative 
angles) of three massless Born momenta is uniquely determined by the T and Tm values. There 
are three essentially different kinematical configurations of the Born system, which we denote 
by 5 (5 = 1,2,3), when P$ is the gluon momentum. In what follows we concentrate?] on the 
most probable configuration with 6 = 3 shown in Fig. [I], in which the gluon belongs to the least 
energetic jet with momentum P 3 . 

The definition of the event plane beyond Born level involves momenta of all final particles, 
and the three hard partons with momentum p a , generally speaking, no longer lie in the plane: 

3 n 

Pa = Pa + q a , K mt = ^ \q ax \ + ^ \k ix \, (1.3) 

a=l i=l 

2 The full answer will be given by the sum over the three configurations weighted by the relative 3-jet cross 
section (unless the gluon jet is tagged). 
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Figure 1: The Born configuration 5 = 3 for T = 0.75 and Tm = 0.48. The thrust T and thrust 
major Tm are along the z- and y-axis respectively. The up, down, left and right hemispheres 
(U, D, L, R respectively) are indicated. 



with q a the recoil momenta of the hard partons. Small x-components of the recoil momenta 
enter the value of the observable K out . The longitudinal (in-plane) components q az , q ay are not 
necessarily small due to hard collinear jet splittings which make the final momentum p a a fraction 
of the initial Born P a . As demonstrated in Q, this rescaling of the longitudinal parton momenta 
gets absorbed into the first hard correction to the emission probability of soft gluons, which is 
then resummed and embodied into the radiator. Bearing this in mind, the recoil momenta can 
be treated as small. 

The standard definition of the event plane as the plane formed by the thrust and thrust-major 
axes leads to the four constraints involving the x- and ^-components of particle momenta: 

qix + 22 kix = ' qi y + y^ J kiy = o , 

R R ^ 

<?2x + Qix ■ &{qi v ) + ^2k ix = 0, q 3x + q lx ■ &(-qiy) + ^2k ix = 0. 

U D 

The definition of the right- (R), up- (U) and down- (D) hemispheres is shown in Fig. [I]. Since all 
partons contribute to K ont and participate in the kinematical relations, to resum and exponentiate 
secondary gluon radiation one should resolve the additive constraints ( |1.3j ) and ( |1.4| ) with the 
help of Mellin and Fourier representations. 

Given the complexity of the event plane conditions (|1.4j ), the analytic derivation of the re- 
sulting distributions turns out to be technically rather involved, especially at the SL accuracy 
which is needed to make quantitative predictions. The essential features of the result, however, 
are rather interesting and transparent from the point of view of the physics involved. In this 
letter we shall present the results and discuss their physical origin and properties. The detailed 
analysis to SL accuracy can be found in 
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We calculate the integrated cross section for fixed T, Tm ~ T and T m smaller than a given 
Kout/Q- To SL accuracy, this cross section can be factorized as 

where da^f 1 / dTdTM is the differential Born 3-jet cross section in the parton configuration 5. The 
factor E5 accounts for the soft radiation emitted by the hard qqg system. The first factor is the 
(non-logarithmic) coefficient function. 

In Section 2 we summarize the results of the PT analysis of the i^ out -distribution obtained 
in [|J . Section 3 is devoted to the leading 1/Q NP corrections to the distribution and the mean. 
The logarithmic enhancement in the NP radiation depends on the aplanarity of each of the hard 
partons (qqg) . These aplanarities, in turn, are the result of the recoil against the PT radiation. 
Since the structure of the recoil depends on the direction of the PT gluon, one finds a peculiar 
colour structure of the NP corrections. We give a simple physical explanation of this structure. 
We conclude in Section 4. 



2 Perturbative i^ out -distribution 

At the DL level, the integrated distribution £ in (|1.5|) is given simply by the Sudakov exponent 
S PT (K ou t) ^ exp(-CV-m 2 -?-), C T = Y j C a = 2C F + N Cl (2.1) 

where C a is the colour factor of parton j^a in the given Born configuration. 

With account of subleading effects the coupling starts to run, the scales of individual parton 
contributions acquire different hard correction factors and become geometry-dependent. The 
recoil of hard partons should be taken into account in the observable K out and in the event plane 
kinematics but, as shown in ||, can be neglected at the PT level in the soft radiation matrix 
elements. The result at SL accuracy reads 

Z PT (K ont ) ~ e-^M K o^). S (a s ]n^-) . (2.2) 

V -ft out / 

The first factor is the exponent of the three two-loop parton radiators R a , 

Ra(K') = C a r (/C t , Q7) , r{K£, QD = - [ Q ^ <2k a ) In 9Z , (2 . 3 ) 
where the running coupling corresponds to the physical (bremsstrahlung; CMW) scheme 0. 

In the configuration 5 = 3 of Fig. ^ we have C\ = C 2 = Cp, C 3 = A^" C) and the geometry- 
dependent PT scales Q PT are 

(Q7f={Q7) 2 = ^e-^, [QD'= {Pl ^fK -^. (2.4) 
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The exponential factors here account for the SL hard corrections due to collinear quark (— 3CV) 
and gluon (— /3 ) splittings. Let us note that the characteristic gluon scale Q 3 is proportional 
to the invariant gluon transverse momentum with respect to the qq pair. Notice that when 
the gluon #3 becomes collinear to the quark (antiquark) #2, the scale Q 3 T decreases, and the 
non-Abelian contribution reduces. 



The function S in ( |2.2| ) is a subleading correction factor. It depends on K out via the SL 
function 

- « - In -9-. (2.5) 

dlnA out 7r A out 

To the first order in r' one has [§] 

5 = 1 - In 2 (2d + C 2 + C 3 ) • r'(K^) + (r /2 ) ; (2d + C 2 + C 3 )* =3 = 3C F + N c . (2.6) 

The fact that the contribution of the R-hemisphere jet #1 is twice larger than that of each of 
the L-hemisphere jets #2,3 has a simple kinematical origin. The differential K ont distribution at 
first order in a s is determined by radiation of a single gluon with a small transverse momentum 
k x off the qqg system. 



dinKmt ^{Eaing + (2C 1 + C 2 + C3)ln2} + OK), (2.7) 

where the first term originates from the DL radiator, and the second one from the derivative of 
the SL factor S in (12. 61). So that the distribution can be written as 



d In K n llt 7T 



in full agreement with the kinematical conditions ( |1.4| ). Indeed, when the secondary gluon is 
emitted in the right hemisphere (see Fig. 0d), all three hard partons experience equal recoils, 



k z > 0, k y > : q lx = -k x = q 2x = -q 3x 
k y < : q lx = -k x = -q 2x = q 3x 



K out = A-\k x \. (2.9) 



On the other hand, for the secondary gluon in the left hemisphere (see Fig. ^|c or d), only one 
hard parton recoils against it: 



k z < 0, k y > : q 2x = -k x ] q lx = q 3x = 
ky < : q 3x = -k x ] q lx = q 2x = 



K out = 2- \k x \. (2.10) 



As we shall see in the next section, similar kinematical effects show up in the structure of the 
NP contribution as well. 

( |2.8| ) is an example of the rich information on the colour structure and on the geometry of 
the underlying event provided by this observable. 
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(c) 
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Figure 2: (a) Hard partons in a generic Born configuration, (b) Soft gluon k (the short 
curly stick) is emitted in the right hemisphere. All three hard partons experience equal out-of- 
plane recoil, see ( |2.9|) . White (shadowed) partons have positive (negative) x-components of the 
momentum, (c) Soft gluon k is emitted in the up-left region. According to ( ^.10[ ), only parton 
recoils (shadowed; has a negative x-components of the momentum), (d) The case of k in the 
down-left region. 



3 NP-effects in the distribution and mean 



Similarly to other jet observables, the leading NP correction to the K out distribution originates 
from radiation of small transverse momentum gluons ("gluers") by the qqg antenna. It is im- 
portant that, at the NP level, the parton recoils q a not only contribute to the observable (( |1.3|) , 
(|1.4[)), as in the PT case, but also affect the radiator itself. In the linear approximation in 1/Q, 



the radiator acquires a NP correction in the form 



r np = u . x *pJ2c a ln^- (3.r 

I Qax I 



a=l 



Here the variable v in (|3.1| ) is Mellin conjugated to K oxA , and A NP (of dimension of mass) is a 
standard NP parameter which can be related to the momentum integral of the QCD coupling 
in the infrared region and parameterizes the 1/Q power correction. Qa P ~ Q are the relevant 
scales which have the same geometric structure as the PT scales ( |2.4| ), and differ from those by 
a finite numerical factor [H. 
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The logarithmic enhancement of the NP contribution (|3.1|) is similar to that for the jet 
Broadening observables (see 0). It originates from the logarithmic integration over the gluer 
angle 9 with respect to the direction of the emitting parton p a , starting from large angles 9 = O (1) 
down to the angle of the hard parton 9 m i n ~ \q ax \/Q- Radiation at smaller angles corresponds to 
collinear splitting and does not contribute to the observable due to real-virtual cancellation. 

It is clear from the structure of the integral involving the Mellin exponent e uK ° ut and the 
exponent of the sum of the PT and NP radiators that the correction term being linear in 
v results in a shift of the argument K out of the PT distribution. To find this shift one has to 
perform the q ax integrals, for a given K ont , taking into consideration the phase space constraints 

0- 

To explain the structure of the answer let us introduce the distributions T> a (\q ax \) of the 
parton recoil momenta. For the leading jet #1 it is given by the derivative of the DL Sudakov 
exponent: 

VM = ^-e- R ^~ 1 \ q=\q lx \. (3.2) 
dmq 

The average value of the logarithm in ( |3.1| ) (the logarithm of the inverse quark angle) evaluated 
under the restriction \qi x \ < K ont becomes 



nNP \ rK out i G m> <oNP 

InfM = ^*& / °± \n ( ^V 1 (q)=\n ( ^ + E 1 (K^), 



(3.3) 

For the partons #2 and #3 in the L-hemisphere the situation is more delicate. As we already 
know from the plane constraints, in order to inhibit the recoil of a left-hemisphere parton j^a we 
have to veto radiation both from j^a and the right-hemisphere parton #1. Therefore, 

V a {q) = J- e-^M-aita- 1 ), a = 2,3, (3.4) 
ding 

and the average logarithms become 

InfM =ln^- + £? B (it^) 

|9»»|<ifottt out 



EJKzti = f K ° Ut ^ e -^("- 1 )-^fe- 1 )+^^ou 1 t)+^^out) , 

q 



(3.5) 



2 O 

In the limit of moderately small K outl so that a s In <C 1, these functions can be approxi- 
mately evaluated to give 

Ei(K^) = - J—r^ ~ ^ -g- + O (1) , 

2 V / ^l«s(Q) K out 

O (3 ' 6) 

2V(^a + Ci)a s (Q) Aout 
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Starting from large values oc the functions E decrease with K out . They vanish as 

E oc l/r'(K~^ t ) in the limit of extremely small K out , such that the Sudakov suppression be- 
comes very strong, r' = ^ In ^> 1. We conclude that the average logarithms of the parton 

recoil angles in (3.1) stay at large constant values O (l/yfa£\, given by the the first terms on 
r.h.s. of (3.6), in a broad region 1 > K out /Q 3> exp(— l/y/o£), and follow In -^—^ for still smaller 



l out 

values of K ont . On the basis of these considerations we can obtain the results for the distribution 
and mean. 



3.1 NP correction to (K, 



out 



The NP correction to the mean value of -Kout is given by the z/-derivative at v = of the NP 
radiator in To evaluate {K ont ) it suffices to take the unrestricted average of \n(Q/\q ax \) by 

setting /Cut m (|3.3| ), ( |3.5|) at their maximum values, i.e. E a (l/Q^ p ). The accurate answer that 
includes first subleading corrections reads || 

3 

(-Kout) = (^out) rl + A*" ■ (Kont)^ (Kont)^ = ^ C a E a [ — ^ ) (1 + O (<*)) . (3.7) 



^Tc a E a ( -p 

0=1 Wa 



NP _ 7T / Cl , C*2 C3 



Analytically, 

(K^f"' ~ — ^ ( + - ~* + - = ) + Oil). (3.8) 

2 v / ^(Q)Vv / cT v/cT+c^ WT+v;) ^ 1 ] 

Since combinations of charges appear in the denominators, we remark that the jets do not 
contribute independently to the mean at the NP level. 



3.2 NP shift of the PT ^-distribution 

The NP shift of the distribution contains two structures, the additive and the non-additive pieces: 

3 

S(K'out) = S PT (K ut - X NP SK), 5K = Y,C a (A« + A'J . (3.9) 

a=l 

The additive contribution is 

NP O 
A a (K out ) = \n^ + i;(l + R')+ lE ~ hi-?-. (3.10) 

The other contribution has the following structure (a more accurate expression can be found 
in i): 

A[(K ont ) ^ E^) , 



T 



c 3 



A' 2 (K mt ) ^ -± E 2 (K^ t ) , (3.11) 
A' 3 (# out ) ~ ^ E 3 (K^ t ) . 
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The additive contribution A has a simple origin. It corresponds to the situation where all 
three emitters have recoil momenta of order A out . Therefore, this contribution to the shift 
follows immediately from the NP radiator (|3.1| ) with \q ax \ ~ K out . Such a situation is typ- 
ical of well-developed parton systems with many secondary partons. In the kinematical re- 
gion a s In 2 jp— 3> 1 the additive term A takes over A', and fully dominates in the region 

"out 

r' = 2s In jP- > 1: 

n "out 

5K=(C 1 + C 2 + C 3 )\n-^- + 0(1), ^ln 2 -^-»l. (3.12) 



On the contrary, for a s In 2 <C 1 there are few secondary partons, and the non-additive 

"out 

contribution A' is dominant. To analyse this situation let us relate the NP correction to the 
integrated spectrum £ to the differential distribution and expand the latter to first order in a s 
in the region a s In <C 1 : 

"out 

£_ S PT^_ A NP^.^_^_ A NP_^.^. R > = CTr >. (3 . 13) 

Putting together the A and A' terms in SK we can then represent the answer as 

R'-6K~R'- C T In + R' 2 • (C^ + C 3 A 3 ) + R' 3 ■ (dEi + C 2 E 2 ) , (3.14) 

A out 

where all the A-functions are evaluated at K~^ t . In the region under consideration we can use 
the approximation (|3.6|), 



K(K£) = E a (Q-') - In = (A out )f - In (3.15) 

-tt-out J^out 

where (A" out )^ P is the NP contribution to the mean associated with parton #a in ( |3.7[ ). Assem- 
bling the logarithmic pieces we arrive at 

Q 



R' ■ 5K ~ R[ ■ C T In 



A n 



+ R' 2 ■ ( <A out )f P + C 2 In + (K out )f) (3.16) 

V -K out / 

+ R's ■ ((Kout)T + (A out )f + C 3 In . 

V A out y 

Now we are ready to discuss the meaning of this expression. 

Making notice that the factor R' a corresponds to the radiation of a perturbative gluon off 
the hard parton #a, we see that the first line describes PT emission in the R hemisphere. In 
this case all three hard partons experience equal recoil, \q ax \ = K out /4: (see Fig. ||b), and the 
logarithmic factors in the NP radiator ( |3.1|) are the same for all three gluers. 
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The second and third lines correspond to PT gluon emission off the parton ^2 and #3. In 
the former case \q2 X \ — K ont /2, while the other two partons remain in the plane (qi x = q% x = 0) 
(see Fig. 0c). Now, NP gluer radiation off #2 brings in a logarithmic correction as above. At 
the same time, the logarithms In Q/\q ax \, a = 1,3 are integrated with the corresponding Sudakov 
distributions T> a producing the contributions identical to those to the mean. 



Substituting (K OVLt )^ P from 



an approximate expression for the NP shift becomes 



5K 



(C2 + C3)\/Ci C2C3 



+ 



Cl 



2 



1 



c 7 



111 



Ct \ \JC\ + C*2 \[C~\ 

Q 



7l 



O(l) 



TV 



2v^(g) 
< 1. 



(3.17) 



The general expression for the NP shift ( |3.9| )- (|3~TT| ) interpolates between the two limiting regimes 
(Km and (KW). 



4 Conclusions 



We have presented the results and discussed the main physical ingredients of the QCD analysis of 
the out-of-event-plane momentum distribution in 3-jet events. It aims at improving the accuracy 
of the theoretical description of the final state structure of hard processes involving more than 
two jets. It should be possible to generalize the present approach to analyse production of two 
large-p_i_ jets in DIS and production of a jet recoiling against W^, Z° or a large-p_i_ photon in 
hadron-hadron collisions. 

The most interesting feature of the K ont distribution is the non-additivity of contributions due 
to the radiation off the three jets. In the PT contribution it occurs at the level of a subleading 
SL correction, while in the 1/Q NP correction terms it shows up at the leading level via peculiar 

combinations of colour charges like (C F + C A )-\ (see @), These can be traced back 



to the specific pattern of parton recoil due to the geometry of the event plane. 

Technical details of the calculations and the final expressions within SL accuracy can be found 
in The relative O (a s ) correction in ( |1.5|) has been computed by matching the expansion 



of the approximate resummed PT formula from M with the O (a^) calculation of the differential 
K ont distribution provided by four-parton event generators. (It is possible, in principle, to further 
improve the result by including the O (a^) term in the non-logarithmic correction factor with 
the help of five-parton generators.) 

Predictions for the K out distribution and means in qq^ events can be obtained from the above 
formulae by simply setting to zero the colour charge corresponding to the gluon (C3 = in the 
colour configuration of Fig. 1.) 

Along the same lines the PT and NP predictions can be derived for the mean and the distri- 
bution of /Cut in the i?-hemisphere, rather than -Kout of the event as a whole. The corresponding 
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perturbative formulae are much simpler, and so are the NP expressions. They can be essentially 
obtained by setting Ci = C3 = in the above expressions for the total K ont . (There is an addi- 
tional difficulty, though, in calculating the proper scales Q FT and Q NP with subleading accuracy, 
due to kinematical complications.) 

Acknowledgements We are grateful to Gavin Salam for discussions and suggestions. 
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